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LONG-TERM  GOALS 

The  long-term  goal  of  this  project  is  to  develop  a  robust  and  hardened  high-resolution  air-ocean 
coupled  tropical  cyclone  (TC)  data  assimilation  and  prediction  system  that  is  able  to  assimilate  the 
wide  variety  of  available  in-si tu  and  remotely-sensed  observations  in  order  to  analyze  and  predict  TC 
structure  and  intensity  changes  in  an  operational  environment.  Significant  gains  have  been  made  in  TC 
track  prediction  over  the  past  three  decades.  This  considerable  achievement  is  due,  in  large  part,  to  the 
steady  improvement  of  numerical  models,  especially  the  global  scale  prediction  systems,  and  the 
judicial  utilization  of  multi-model  ensemble  results.  In  contrast,  the  TC  intensity  forecast  by  numerical 
models  has  shown  very  little  improvement  during  the  same  time  period,  and  remains  a  formidable 
forecast  problem.  Advanced  statistical  prediction  models  nowadays  are  able  to  predict  the  trend  for 
intensification,  but  as  statistical  tools,  they  inherently  cannot  predict  the  rapid  intensity  changes,  as 
evident  in  Katrina  of  2005,  and  Roke  in  the  W.  Pacific  in  20 1 1 .  It  is  generally  accepted  now  that  while 
advancements  in  data  assimilation  and  modeling  have  resulted  in  better  analyses  and  predictions  of 
steering  flows,  the  processes  that  affect  the  structure  and  intensity  of  tropical  cyclones  are  much  more 
difficult  for  current  numerical  models  to  capture  and  reproduce.  Physical  processes  in  tropical 
cyclones  that  can  affect  their  structure  and  intensity  include  enthalpy  and  mechanical  interchanges  with 
the  underlying  ocean  and  land  surfaces,  shallow  and  deep  atmospheric  convection  in  the  convectively 
unstable  tropical  atmosphere  with  vertical  and  horizontal  wind  shears,  and  internal  multi-scale  non¬ 
linear  dynamic  interactions.  Current  prediction  systems  have  been  shown  to  be  able  to  reproduce  rapid 
intensification  in  case  studies  involving  complex  upper  tropospheric  and  oceanic  conditions  in  a 
carefully  conducted  simulation  model. 

OBJECTIVE 

The  objective  is  to  develop,  validate,  and  transition  new  methods  to  initialize  a  coupled  tropical 
cyclone  (TC)  modeling  system  using  advanced  data  assimilation  and  advanced  physical 
parameterization  methods  in  order  to  achieve  superior  short  term  forecasts  (0-2  days)  of  TC  structure 
and  intensity.  Our  approach  is  to  leverage  emerging  data  assimilation,  modeling,  and  air-ocean 
coupling  techniques,  as  well  as  observational  results  from  the  scientific  community  including  the  recent 
T-PARC/TCS08  and  ITOP  field  campaigns  to  build  upon  the  existing  modeling  capabilities.  We  will 
develop  multiple  approaches  to  initialize  the  tropical  cyclone  including  (from  least  to  most 
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sophisticated  and  challenging  method):  i)  bogus  initialization,  ii)  three-dimensional  variational  analysis 
(3D-Var)  using  synthetic  observations,  iii)  dynamical  initialization  (TCDI),  iv)  an  ensemble  Kalman 
filter  approach,  and  v)  four-dimensional  variational  (4D-Var)  analysis.  The  boundary  layer  and  cloud 
microphysics  parameterizations,  which  have  a  direct  impact  on  the  short-term  forecasts  as  well  as  the 
analysis  background  fields,  will  be  developed  further  and  new  methods  will  be  tested. 

APPROACH 

Our  approach  is  to  integrate  emerging  data  assimilation,  modeling,  and  coupling  techniques,  as  well  as 
observational  results  from  recent  Office  of  Naval  Research  (ONR)  field  campaigns,  into  the  existing 
framework  in  the  Coupled  Ocean/ Atmosphere  Mesoscale  Prediction  System  (COAMPS®1)  for 
applications  to  the  analysis  and  prediction  of  TC  position,  structure,  and  intensity  (referred  to  as 
COAMPS-TC).  Our  specific  goal  is  to  improve  the  0-48  h  COAMPS-TC  analysis  and  prediction  of 
intensity  and  structure.  Specific  technologies  that  will  be  developed,  tested,  and  integrated  into 
COAMPS  for  application  to  TC  prediction  in  this  project  using  COAMPS-TC  are:  (1)  TC  initialization 
and  analysis  techniques,  (2)  TC  physical  process  parameterizations  with  particular  emphasis  on  the 
microphysics,  and  (3)  air-ocean-wave  coupling  techniques.  The  proposed  project  will  also  leverage 
other  research  that  has  taken  place  on  COAMPS-TC  including  the  recently  completed  COAMPS-TC 
RTP  and  the  coupled  modeling  development  through  the  Battlespace  Environments  Institute  (BEI), 
which  addressed  the  coupling  of  the  COAMPS  atmospheric  model  to  the  Navy  Coastal  Ocean  Model 
(NCOM)  and  the  SWAN  and  WWIII  wave  models  using  the  Earth  System  Modeling  Framework 
(ESMF).  Also,  the  project  will  leverage  emerging  scientific  findings  from  the  recent  TCS08  and  ITOP 
in  the  Northwestern  Pacific. 

WORK  COMPLETED 

We  are  currently  developing  and  evaluating  several  approaches  to  initialize  the  tropical  cyclone  vortex 
and  environment  including  (from  least  to  most  sophisticated  and  challenging  method):  i)  bogus 
initialization,  ii)  three-dimensional  variational  (3D-Var)  analysis  using  synthetic  observations,  iii) 
dynamical  initialization  (TCDI),  iv)  an  ensemble  Kalman  filter  (EnKF)  approach,  and  v)  four¬ 
dimensional  variational  (4D-Var)  analysis.  Versions  of  several  these  methods,  namely  the  3D  Var  with 
synthetics  and  TCDI,  have  already  been  developed  for  COAMPS-TC  and  to  varying  degrees  evaluated. 
However,  their  performance  has  been  lacking  in  a  number  of  aspects  including  short  term  intensity  and 
track  error  issues.  Greater  promise  for  more  accurate  initialization  of  the  TC  vortex  exists  for  the  more 
sophisticated  EnKF  and  4D-Var  methods.  The  EnKF  and  4D-Var  methods  will  be  advanced  further  as 
part  of  this  project,  and  all  of  the  methods  will  be  statistically  evaluated  for  a  large  set  of  TC  cases 
using  a  variety  of  datasets  include  remote  sensing. 

COAMPS-TC  has  been  tested  in  both  coupled  and  uncoupled  modes  over  the  past  several  tropical 
cyclone  seasons  in  the  Pacific  and  Atlantic  basins.  In  these  applications,  the  atmospheric  portion  of  the 
COAMPS-TC  system  makes  use  of  horizontally  nested  grids  with  grid  spacing  of  45,  15,  and  5  km. 

The  15-  and  5-km  grid-spacing  meshes  track  the  TC  center,  which  enables  the  TC  convection  to  be 
more  realistically  represented  on  the  finest  mesh  in  an  efficient  manner.  The  forecasts  make  use  of  the 
Navy  or  the  NOAA  global  model  for  lateral  boundary  conditions.  The  forecasts  are  routinely 
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disseminated  in  real  time  to  NHC,  JTWC,  and  HFIP  researchers.  The  forecast  graphics  are  also 
available  in  real  time  at  http://www.nrlmry.navy.mil/coamps-web/web/tc. 

A  two-stage  dynamic  initialization  scheme  for  COAMPS-TC  has  been  developed  and  evaluated.  The 
scheme  consists  of  a  tropical  cyclone  dynamic  initialization  (TCDI)  method  in  conjunction  with  a 
forward  dynamic  initialization  (DI).  In  the  TCDI  stage,  a  dynamically  spun-up  vortex  matching  the 
observed  intensity  is  inserted  into  the  COAMPS-TC  initial  conditions  after  3D-Var.  In  the  subsequent 
DI  stage,  COAMPS-TC  is  integrated  for  a  period  of  12  h,  with  Newtonian  relaxation  to  the  modified 
3D-Var  analysis  horizontal  momentum  with  the  addition  of  the  TCDI  vortex.  This  two  stage  dynamic 
initialization  scheme  improves  both  the  vortex  nonlinear  balance  and  consistency  in  the  transverse 
circulation,  while  reducing  spurious  gravity  wave  activity  generated  by  interpolation  and  analysis  errors 
in  the  3D-Var  scheme.  In  addition  the  benefits  listed  above,  this  two-stage  dynamic  initialization 
scheme  represents  only  a  minor  additional  cost  to  the  3D-Var  scheme,  and  it  is  significantly  cheaper  in 
a  computational  sense  than  either  the  4D-Var  or  EnKF  methods. 

RESULTS 

a.  COAMPS-TC  Real-Time  and  Retrospective  Model  Forecasts 

Real-time  COAMPS-TC  forecasts  have  been  conducted  using  U.S.  Department  of  Defense  High 
Perfonnance  Computing  (HPC)  platforms  over  the  past  several  years.  An  example  of  the  intensity 
forecast  perfonnance  of  COAMPS-TC  for  a  large  number  of  cases  (more  than  450  cases  at  the  24  h 
forecast  time)  in  the  W.  Atlantic  region  for  the  2010  and  2011  seasons  is  shown  in  Figure  1  (for  a 
homogeneous  statistical  sample).  The  COAMPS-TC  model  had  the  lowest  intensity  error  of  any 
dynamical  model  for  the  36-120  h  forecast  times,  which  is  an  important  period  for  forecasters  and 
decision  makers.  Other  numerical  models  included  in  this  analysis  are  operational  models  run  by 
NOAA  (HWRF,  GFDL),  and  the  Navy’s  current  operational  limited  area  model  (GFDN).  This 
promising  performance  is  a  result  of  a  large  effort  devoted  to  developing  and  improving  COAMPS-TC 
over  the  past  several  years.  Key  aspects  of  the  COAMPS-TC  system  that  have  contributed  most  to  the 
forecast  intensity  skill  improvement  (as  deduced  from  sensitivity  tests)  include:  i)  new  synthetic 
observations  and  data  assimilation  system  modifications;  ii)  modified  turbulence  parameterization  (in 
particular  the  mixing  length)  in  and  above  the  boundary  layer;  iii)  surface  drag,  heat,  and  moisture 
exchange  consistent  with  recent  field  campaigns;  and  iv)  fidelity  of  the  microphysical  parameterization. 

An  example  of  a  real-time  COAMPS-TC  forecast  for  the  recent  Hurricane  Irene  is  shown  in  Figure  2. 
The  composite  National  Weather  Service  radar  reflectivity  is  shown  in  the  top  panel  near  the  time  of 
landfall  in  North  Carolina  at  1 148  UTC  27  August  2011  and  the  COAMPS-TC  predicted  radar 
reflectivity  at  36  h  valid  at  1200  UTC  is  shown  in  the  bottom  panel.  The  COAMPS-TC  forecast  shown 
in  Figure  2  is  for  the  model  second  grid  mesh  (15  km  horizontal  resolution).  The  model  prediction  was 
accurate  in  the  track  (skill  similar  to  other  dynamical  models),  eventual  landfall  location  (Figure  2), 
storm  intensity  (Figure  4),  as  well  as  the  structure  and  size  (from  a  qualitative  perspective),  an 
especially  important  characteristic  of  this  particular  storm  in  such  close  proximity  to  the  U.S.  East 
Coast.  One  noteworthy  aspect  of  Irene  was  its  large  size,  with  tropical  stonn  force  winds  (34  knot) 
extending  radially  outward  from  the  eye  for  nearly  300  km.  The  large  size  of  Irene  is  also  apparent  in 
the  observed  radar  reflectivity  in  Figure  2.  The  COAMPS-TC  prediction  captures  the  large  areal  extent 
of  the  precipitation  field,  as  well  as  its  asymmetry  about  the  TC  center  (most  of  the  precipitation  is 
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north  and  east  of  the  center.  This  large  shield  of  heavy  precipitation  caused  severe  river  flooding  as  it 
slowly  moved  north  through  the  mid-Atlantic  and  Northeast  U.S.  The  simulated  radar  reflectivity  for 
the  COAMPS-TC  grid  mesh  3  (5  km  horizontal  resolution),  shown  in  Figure  3,  illustrates  the 
capability  of  the  model  to  capture  the  (incr-scale  features,  such  as  the  eyewall  and  rainbands,  in 
generally  good  agreement  with  the  observed  reflectivity.  Sensitivity  tests  (not  shown)  indicate  that  the 
synthetic  observations  and  TC-related  data  assimilation  system  modifications,  along  with  the  in-cloud 
turbulent  mixing  representation,  are  important  for  the  proper  simulation  Irene’s  structure.  It  should 
also  be  noted  that  COAMPS-TC  makes  use  of  a  1.5  order  closure  turbulence  parameterization  that 
predicts  turbulent  kinetic  energy  (TKE)  and  includes  advection  of  TKE  as  well.  We  hypothesize  that 
the  advection  of  TKE  may  be  important  for  the  proper  representation  of  the  turbulence  evolution  within 
the  hurricane  boundary  layer  beneath  and  near  the  eyewall. 

Overall,  the  Navy's  COAMPS-TC  real-time  intensity  predictions  of  Hurricane  Irene  were  more  skillful 
than  the  other  leading  operational  governmental  forecast  models,  as  shown  in  Figure  4.  All  of  the 
available  models  except  for  COAMPS-TC  had  a  tendency  to  over-intensify  Irene,  often  by  a  full  stonn 
category  or  more.  These  real-time  COAMPS-TC  forecasts  were  used  by  forecasters  at  the  National 
Hurricane  Center  as  part  of  an  experimental  HFIP  multi-model  ensemble.  The  COAMPS-TC 
consistently  provided  accurate  real-time  intensity  forecasts  during  the  period  23-28  August  2011,  when 
critical  decisions  were  made  by  forecasters  and  emergency  managers  including  evacuations. 

b.  COAMPS-TC  Ensemble  Forecasts  Experiments 

A  new  COAMPS-TC  ensemble  system  that  is  capable  of  providing  probabilistic  forecasts  of  TC  track, 
intensity,  and  structure  has  been  developed  by  scientists  at  NRL  in  Monterey,  CA.  This  system  makes 
use  of  the  community-based  Data  Assimilation  Research  Testbed  (DART)  capability  developed  at  the 
National  Center  for  Atmospheric  Research,  which  includes  various  options  for  Ensemble  Kalman  Filter 
(EnKF)  data  assimilation.  The  COAMPS-TC  DART  system  constitutes  a  next  generation  data 
assimilation  system  for  tropical  cyclones  that  uses  flow  dependent  statistics  from  the  ensemble  to 
assimilate  observational  information  on  the  mesoscale. 

A  real-time  COAMPS-TC  ensemble  system  was  run  in  a  demonstration  mode  in  201 1  for  the  W. 
Atlantic,  E.  Pacific,  and  W.  Pacific  regions.  The  system  comprised  an  80-member  COAMPS-TC 
cycling  data  assimilation  ensemble  on  three  nested  grids  with  horizontal  spacing  of  45,  15,  and  5  km. 
Ten-member  forecasts  were  perfonned  twice  daily  to  five  days  using  the  same  three  nested  grid 
configuration  as  the  data  assimilation  ensemble.  The  first  10  members  of  the  data  assimilation 
ensemble  were  used  to  define  the  forecast  ensemble.  Lateral  boundary  conditions  for  the  forecast 
ensemble  were  drawn  from  the  GFS-EnKF  ensemble  forecast.  Examples  of  probabilistic  products  for 
Hurricane  Irene  are  shown  in  Figure  5  for  both  track  (top  panel)  and  intensity  (bottom  panel).  This  is  a 
real  time  forecast  initialized  at  1200  UTC  23  August,  which  is  four  days  prior  to  landfall.  The 
probabilistic  track  product  shows  the  TC  position  from  the  individual  ensemble  members  every  24  h 
and  ellipses  that  encompass  1/3  and  2/3  of  the  ensemble  member  forecast  positions.  Note  that  the 
observed  landfall  location  of  the  eye  (see  Figure  2)  was  within  the  ensemble  distribution,  although  the 
ensemble  mean  landfall  was  approximately  12  h  later  than  observed.  The  probabilistic  intensity 
product  (lower  panel)  shows  a  considerable  spread  among  the  members,  particularly  beyond  84  h,  just 
prior  to  landfall.  These  products  can  be  extremely  valuable  to  assess  the  uncertainty  in  both  track  and 
intensity  forecasts,  and  NRL  is  currently  developing  these  capabilities  and  products  further. 
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c.  Coupled  COAMPS-TC  Forecast  Experiments 

The  COAMPS-TC  system  was  run  in  a  fully  coupled  mode,  interactive  with  NCOM,  during  the  Office 
of  Naval  Research  sponsored  Interaction  of  Typhoon  and  Ocean  Project  (ITOP)  during  the  summer  and 
fall  of  2010.  An  example  of  a  fully  coupled  COAMPS-TC  forecast  for  Typhoon  Fanapi  perfonned  in 
real  time  is  shown  in  Figure  6.  The  NCOM  ocean  model  was  applied  using  a  10-km  horizontal  grid 
increment  in  this  example,  and  the  finest  mesh  for  the  atmospheric  model  used  a  5 -km  grid  increment. 
The  atmosphere  and  ocean  fluxes  were  exchanged  every  ocean  model  time  step.  The  COAMPS-TC 
predicted  track  (red)  from  a  90-h  real  time  forecast  valid  at  0600  UTC  19  September  2010  is  quite 
close  to  the  observed  or  best  track  (black).  The  sea  surface  temperature,  shown  in  color  shading, 
indicates  significant  cooling  was  predicted  by  COAMPS-TC  during  the  passage  of  Fanapi  due  to 
enhanced  mixing  by  the  strong  near-surface  winds  of  the  typhoon.  The  predicted  cooling  of  the  sea 
surface  temperatures  of  2-4°C  is  in  agreement  with  estimates  from  in  situ  and  remote  sensing 
observations  in  this  region.  The  negative  impact  of  the  ocean  cooling  underneath  the  tropical  cyclone 
can  reduce  the  TC  intensity  and  broaden  the  tropical  cyclone  secondary  circulation,  which  underscores 
the  importance  of  properly  representing  these  air-sea  interaction  process. 

d.  Dynamic  Initialization  Development  and  Experiments 

An  evaluation  of  the  track  and  intensity  perfonnance  of  the  two-stage  dynamic  initialization  scheme 
was  conducted.  Historical  TCs  were  simulated  using  COAMPS-TC  with  the  3D-Var  and  synthetics 
initialization  method  (CNTL)  and  the  two-stage  dynamic  initialization  method  (TCDPDI). 
Improvements  in  both  the  initial  vortex  dynamic  and  thermodynamic  balance  and  intensity  and  track 
errors  were  found  by  using  TCDI/DI.  An  example  of  an  improvement  in  the  initial  balance  is  shown  in 
Fig.  7.  As  shown  by  Fig.  7,  in  the  CNTL  initialization  method,  a  spin-down  occurs  early  in  the  forecast 
causing  the  large  errors  in  the  subsequent  intensity.  In  the  TCDI/DI  method,  the  TC  is  in  a  state  of 
improved  initial  balance,  leading  to  more  minor  adjustments  early  and  a  much  improved  intensity 
forecast  in  the  first  6  hours.  Both  initializations  were  tested  on  multiple  TCs  in  the  western  North 
Pacific  and  North  Atlantic  basins  in  2010-2011.  The  intensity  and  track  errors  on  a  large  sample  of 
cases  are  given  in  Fig.  8.  As  shown  by  Fig.  8,  the  TCDI/DI  method  reduced  the  average  track  errors 
slightly  (approximately  10-30  nm,  left  panel)  and  reduced  the  average  intensity  errors  significantly 
(approximately  25%,  with  statistical  significance  at  many  lead  times,  right  panel).  These  results 
demonstrate  that  the  TCDFDI  method  can  be  used  to  improve  TC  intensity  forecasts.  This  method 
combines  the  benefits  from  the  TCDI  method  (red  curve)  the  early  part  of  the  forecast,  and  the  DI 
method  (purple)  curve  later  in  the  forecast  to  produce  improvements  at  all  lead  times,  in  comparison  to 
the  CNTL  method. 

TRANSITIONS 

The  tropical  cyclone  application  of  COAMPS  will  transition  to  6.4  projects  within  PE  0603207N 
(PMW-120)  that  focus  on  the  transition  COAMPS  to  FNMOC. 

RELATED  PROJECTS 

COAMPS  will  be  used  in  related  6. 1  projects  within  PE  0601 153N  that  include  studies  of  air-ocean 
coupling  and  boundary  layer  studies,  and  in  related  6.2  projects  within  PE  0602435N  that  focus  on  the 
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development  of  the  atmospheric  components  (QC,  analysis,  initialization,  and  forecast  model)  of 
COAMPS. 
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Figure  1.  Wind  speed  mean  absolute  error  (MAE)  (knots;  1  knot=0.514  m  s'1)  as  a  function  of 
forecast  time  for  the  2010  and  2011  seasons  in  the  Atlantic  basin  for  a  homogeneous  statistical 
sample.  The  numerical  models  included  in  this  analysis  are  the  Navy’s  COAMPS-TC,  operational 
models  run  by  NOAA  (HWRF,  GFDL),  and  the  Navy’s  current  operational  limited  area  model 
(GFDN).  The  number  of  cases  is  shown  at  the  bottom. 
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Figure  2.  The  NWS  composite  radar  reflectivity  (from  NOAA)  valid  at  1148  UTC  27  August  2011 
(left)  (source  NOAA)  and  the  COAMPS-TC  3  6-h  forecast  radar  reflectivity  performed  in  real  time 
and  valid  at  1200  UTC  27  August  (right)  for  Hurricane  Irene.  The  COAMPS-TC  reflectivity  is 
shown  for  the  second  grid  mesh,  which  has  a  horizontal  resolution  of  15  km. 
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Figure  3.  The  COAMPS-TC  3  6-h  forecast  radar  reflectivity  valid  at  1200  UTC  27  August  for  the 
third  grid  mesh,  which  has  a  horizontal  resolution  of  5  km. 


Intensity  error,  NHC  criteria 


Figure  4.  Wind  speed  mean  absolute  error  (MAE)  (knots)  as  a  function  of forecast  time  for 
Hurricane  Irene  for  a  homogeneous  statistical  sample.  The  numerical  models  included  in  this 
analysis  are  the  Navy’s  COAMPS-TC,  operational  models  run  by  NOAA  (HWRF,  GFDL),  and  the 
Navy’s  current  operational  limited  area  model  (GFDN).  The  number  of  cases  is  shown  at  the 
bottom.  Only  forecasts  after  Irene  has  moved  away  from  Hispaniola  are  shown  here. 
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TC  =  OBL,  DTG  =  2011082312 


Figure  5.  Probabilistic  products  from  the  COAMPS-TC  ensemble  for  Hurricane  Irene 
corresponding  to  the  track  (top  panel)  and  hurricane  intensity  (bottom  panel).  This  real  time 
forecast  was  initialized  at  1200  UTC  23  August,  which  is  approximately  four  days  prior  to  landfall. 
The  probabilistic  track  product  shows  the  TC position  from  the  individual  ensemble  members  every 
24  h  and  ellipses  that  encompass  1/3  and  2/3  of  the  ensemble  members.  The  intensity  (knots) 
distribution  is  shown  as  a  function  of  forecast  lead  time  (hours)  with  the  minimum  value,  maximum 
value  and  various  quantiles  of  the  ensemble  distribution  shown  as  denoted  by  the  legend. 
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Figure  6.  The  COAMPS-TC predicted  track  (red)  for  Typhoon  Fanapi  from  a  90  h  real-time 
forecast  valid  at  0600  UTC 19  September  2010  and  the  observed  best  track  (black).  The  sea  surface 
temperature  (°C)  is  shown  in  color  shading  and  indicates  significant  cooling  during  the  passage  of 
Fanapi.  The  surface  currents  are  shown  by  the  white  vectors. 
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Figure  7.  The  intensity  max  maximum  sustained  wind  (VMAX,  m/s)  of  both  the  CNTL  and 
TCDI/DI  runs  in  the  first  6  hours  of  the  COAMPS-TC  forecast  of  16W  (Bolaven)  in  the  western 
North  Pacific  basin,  in  comparison  to  the  best  track  estimate  of  the  intensity  from  the  Joint  Typhoon 
Warning  Center  (JTWC).  Saffir-Simpson  hurricane  categories  are  also  given. 
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2010-2011  Large  Sample  TRACK  ERROR  2010-2011  Large  Sample  INTENSITY  ERROR 


Figure  8.  (left  panel)  The  average  track  errors,  and  (right  panel)  the  average  intensity  errors  for  a 
large  sample  of  cases  (120  at  the  initial  time  to  59  at  the  120  h  lead  time).  The  CNTL  is  depicted  by 
the  black  curve,  and  TCDI/DI  is  depicted  by  the  blue  curve.  Additional  curves  for  running  DI 
(purple)  and  TCDI  (red)  alone  are  also  given.  In  the  left  panel,  the  total  track  error  is  solid,  the 
along-track  error  is  dotted,  and  the  cross-track  error  is  dashed.  In  the  right  panel,  the  average 
errors  are  solid,  and  the  biases  are  dotted.  Black  diamonds  indicate  that  a  statistically  significant 
difference  exists  between  any  of  the  colored  error  curves  and  the  black  CNTL. 
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